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Abstract 

We investigate properties of non-translation-invariant measures, describing parti- 
cle systems on Z, which are asymptotic to different translation invariant measures 
on the left and on the right. Often the structure of the transition region can only 
be observed from a point of view which is random — in particular, configuration 
dependent. Two such measures will be called shift equivalent if they differ only 
by the choice of such a viewpoint. We introduce certain quantities, called trans- 
lation sums, which, under some auxiliary conditions, characterize the equivalence 
classes. Our prime example is the asymmetric simple exclusion process, for which 
the measures in question describe the microscopic structure of shocks. In this 
case we compute explicitly the translation sums and find that shocks generated 
in different ways — in particular, via initial conditions in an infinite system or by 
boundary conditions in a finite system — are described by shift equivalent mea- 
sures. We show also that when the shock in the infinite system is observed from 
the location of a second class particle, treating this particle either as a first class 
particle or as an empty site leads to shift equivalent shock measures. 
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1. Introduction 

A major component of statistical mechanics, especially its mathematical aspect, is the 
study of measures or probability distributions for infinite particle systems. Such infinite 
systems represent idealizations of macroscopic physical systems whose spatial extension, 
although finite, is very large on the microscopic scale of interparticle distances or inter- 
actions. The advantage of this idealization is that many phenomena which are clearly 
manifested in real macroscopic systems, such as phase transitions, have precise counter- 
parts in the behavior of the infinite volume measures. The inevitable boundary and finite 
size effects present in real systems, which are frequently irrelevant to the phenomena of 
interest, are eliminated in the thermodynamic (infinite volume) limit [1]. 

Our mathematical characterization of these measures for infinite particle systems is 
very good for situations in which the measures are translation invariant (TI) [1,2]. The sit- 
uation becomes less transparent when dealing with spatially nonuniform measures. These 
can arise in various ways. A rather trivial case occurs when the interaction Hamiltonian or 
the dynamics specifying the evolution is position dependent. More interesting cases arise 
when the translation symmetry is broken "spontaneously" by the measure. We first illus- 
trate these by a well known example from equilibrium statistical mechanics, then discuss 
a nonequilibrium example which is the main focus of this paper. 

Consider the Gibbs measures for the nearest neighbor ferromagnetic Ising model on 
Z d at low temperatures. In d > 3 there exist, in addition to the TI extremal measures fi + 
and //_, in which the spontaneous magnetizations are ±m* with m* ^ 0, many non-TI 
measures called Dobrushin states; a family of these can be obtained as the infinite volume 
limit of systems with ± boundary conditions in the ei-direction [3]. Specifically, let the 
domain O containing the system consist of sites j = {ji, ■ ■ • ,jd) such that j± G Z and 
32, ■ ■ • , 3d £ [-N, -/V] d_1 , with all spins outside O being equal to +1 for ji > and —1 for 
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ji < 0. In each system configuration we consider the set of (d — 1) -dimensional surfaces 
separating +1 and —1 spins, formed of (d — 1) -dimensional faces of cubes in the dual 
lattice 7h d + (1/2, . . . , 1/2). The Dobrushin interface is the maximal connected component 
of this set which contains all such (d — l)-faces outside O. If d > 3 then at low enough 
temperatures (below the roughening transition) this interface remains localized near the 
ji = plane as iV increases. Consequently, in the resulting infinite volume states the 
expectation values (or correlations) depend on the e\ coordinate, e.g., if aj = ±1 is the 
spin at site j G Z d then (o"j) M (±) is positive for j\ > and negative for j\ < [3]. This 
non-translation-invariant infinite volume Gibbs state is one of an infinite family obtained 
via translations in the e\ direction. 

If d = 2, the same boundary conditions produce a translation invariant state in the 
infinite volume limit. This is because the interface, while remaining locally sharp, fluctuates 
in position with the variance of its displacement from the plane j\ = growing like N [4] . 
Consequently, the limiting measure (defined by the N — > oo limit of local correlation 
functions) is a superposition, with equal weights, of the extremal translation invariant 
measures fi+ and //_ [4,5]. Suppose, however, that we view the system from some point 
attached to the Dobrushin interface; for example, we might choose the point (j*, 0), where 
j* is as large as possible so that (j* — 1/2, 0) is on the interface. Note that the value of j* 
and hence the viewpoint will depend on the configuration under view. It seems clear that 
when N — > oo a limiting measure will exist which will not be translation invariant, but 
will instead approach the state (respectively as one goes to infinity in the positive 
(respectively negative) e\ direction; this has not been explicitly established but for results 
in this direction see [6] . Other viewpoints are of course possible, and the resulting measure 
will depend in a complicated way on the choice made. One might choose, for example, 
(j* + ji 1 32) for some fixed (j\ , j'2) , or (j* , 0) with j* as small as possible so that (j* + 1/2, 0) 
is on the interface. One could even add an additional randomness by choosing either (j*, 0) 
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or (j*,0) with equal probability; this seems artificial in the current context but this sort 
of additional randomness is natural and necessary in the one-dimensional system to be 
studied shortly. 

This example illustrates two ways in which non-TI measures arise. The non-TI mea- 
sures for d > 3 arise in viewing the system from nonrandom, fixed, frames. Choice of a 
different frame simply effects a translation of the measure. In contrast, the non-TI mea- 
sures for d = 2 can be seen only if one views the system from a random position — random 
in the sense that it depends on the configuration. Moreover, since the choice of a view- 
point is rather arbitrary and since the effect on the measure of a change in viewpoint is 
hardly transparent, one must now consider a large family of distinct measures arising from 
different viewpoints. Of course, even in d > 3 we could consider the measure as seen from 
a point attached to the interface. At low temperature there seems to be little to be gained 
from such an approach, but it might be of interest between the roughening and critical 
temperatures in d = 3, where the situation is expected to be similar to that in d = 2. 

Several questions arise when the transition region must be described by non-TI mea- 
sures obtained from configuration-dependent viewpoints. Is there a natural choice of view- 
point which will give a best or simplest description of the local structure of the transition 
region? How can one extract intrinsic properties of this region from such a description- 
properties independent of the choice of viewpoint? And, given two non-TI measures, how 
may one decide if they in fact describe the same system seen from different points of view? 
The purpose of the present work is to address such questions. Our motivation, and the 
focus of our study, is in fact not the above example but a nonequilibrium system, the one 
dimensional asymmetric simple exclusion process (ASEP). The results are, for the moment, 
also specific to one-dimensional systems. 

We now describe the non-TI measures arising in the ASEP. The latter [7,8] is a 
model of particles moving on the lattice Z; a configuration r\ of the system has the form 
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rj = (r}(i))n=.z, where rj(i) is or 1 at an empty or occupied site, respectively. Dynamically, 
each particle attempts to jump to a neighboring site, at random times with rate 1, choosing 
its right or left neighbor with probabilities p and q, respectively, where p > 1/2 and 
q = 1 — p. The jump takes place if and only if the target site is empty. The extremal 
stationary TI states of this system are the product (Bernoulli) measures v p with constant 
density p satisfying < p < 1 [8]. There exist also non-TI stationary states, which are 
product measures with nonuniform density approaching as i — > — oo and 1 as i — > oo. 
The latter are in fact special examples of the general class of non-TI measures we will study 
here: those which describe the microscopic structure of shocks present in the macroscopic 
description of the ASEP. 

The ASEP is described on the (Euler) macroscopic space-time scale by the inviscid 
Burgers equation for the particle density t) G [0, 1], where ([9-11]): 

^ + (p-q)-^n(l-n) = 0. (1.1) 
Equation (1.1) has shock solutions, n{x,t) = u{x — Vi), where 

»<»> = {£ till w 

here p+ > p- and the velocity is V = (p — q)(l — p + — p_). A natural question then 
is what behavior of the ASEP system on the microscopic level corresponds to this shock 
solution. For example, one may take the initial state po of the system to be a product 
measure with density at site j given by p- for j < and p + for j > 0, and ask about the 
t — > oo limiting behavior of the state p t at time t. It might seem that if one were to view 
the system from a frame moving with the shock velocity V then one would see in this limit 
a non-TI state describing the intrinsic microscopic structure of the shock. But this is not 
true: because fluctuations in the shock position become unbounded on the microscopic 
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scale as t — > oo, the resulting measure is an equal superposition of the product measures 
u p+ and v p _ [12,13]. It has been shown [14,15], however, that there exists a (nonunique) 
time-dependent random position X t such that the t — > oo limit of the measure fit seen 
from the viewpoint X t , which we shall denote by fx', exists and is spatially asymptotic to 
the product measures u p+ and v p _ : 

lim T-V = i/ P± . (1.3) 

Here T is the translation operator, which acts on configurations by {Trf)(i) = n(i — 1), 
on functions of configurations by (Tf)(rj) = /(T -1 ??), and on measures on configurations 
space by (f)r^ = /) jU . The situation is thus analogous to that of the two dimensional 
Ising model: in the t — > oo limit here, and in the N — > oo limit there, one must look from 
a configuration-dependent viewpoint to see the non-TI state. 

The random position X t discussed above is given by the location of a single second 
class particle inserted into the system, which is then treated as an empty site in obtaining 
the measure //. This viewpoint is doubly random, in that the random configuration rj 
does not completely determine the viewpoint, but only its distribution. The measure fi' 
is invariant under the ASEP dynamics for the system seen from the second class particle 
(we will describe this dynamics below). In previous works [16,17] explicit formulas were 
obtained for a measure ju invariant for this same dynamics and with the same spatial 
asymptotics (1.3), and it is this measure that will be our main example here; presumably 
fi, = fj,', although this has not been established. (It is n' which has shown to be obtained 
by the long time asymptotics described above.) 

In this paper we will focus on questions like those raised above in the context of the 
Ising model, which arise from the possibility of different choices of viewpoint. In Section 2 
we describe the evolution of the ASEP with a second class particle and the resulting 
viewpoint on the shock, as well as several other possible choices of viewpoint. In Section 3 
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we formalize, in a general one-dimensional context, the relation of shift equivalence on non- 
TI measures under which equivalent measures differ by a random change of viewpoint. 
There we define also certain quantities, called translation sums, which characterize this 
equivalence: two measures (which must satisfy certain additional conditions) are shift 
equivalent if and only if all the translation sums for the two measures agree. This result 
will be established in a separate paper [18]. 

In the remainder of the paper we apply these general ideas to the ASEP. In Section 4 
we utilize the results of [17] to compute the translation sums explicitly. From this com- 
putation (and using the verification, here omitted, that the ASEP shock measures satisfy 
the additional conditions mentioned above) we establish both negative and positive results 
about the ASEP shock. In Section 5 we show that for certain values of the parameters the 
shock measure is not shift equivalent to any product measure with a monotone density; 
we show also that when the shock is observed from the location of a second class particle, 
treating this particle either as a first class particle or as an empty site leads to shift equiv- 
alent shock measures. Finally, in Section 6 we show that certain shocks arising in versions 
of the ASEP with different boundary conditions are in fact shift equivalent. 
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2. Points of view for the ASEP shock 

In this section we illustrate the nature of shock measures by considering various view- 
points on the shock for the ASEP, beginning with the viewpoint from a second class 
particle. We first describe briefly the properties of the ASEP when a single second class 
particle is introduced into the system [12]. The second class particle has its own dynamics: 
it attempts to jump exactly as does an ordinary (first class) particle, succeeding only if 
the target site is empty; on the other hand, when a first class particle attempts to jump 
onto the site occupied by the second class particle, the jump succeeds and the two particles 
exchange sites. A configuration of this system is r = (r(i))i e z, where r(i) is if site i 
is unoccupied, 1 if it is occupied by one of the original particles, now called first class 
particles, and 2 if it is occupied by the second class particle. Let us denote the location 
of the second class particle by X. If At is a measure on this system evolving under the 
above dynamics and X t the corresponding location of the second class particle at time t, 
we write T~ Xt X t for the measure describing the configurations as seen from the second 
class particle. 

There is an alternate, equivalent way to describe the system with a second class particle 
[12]. Consider two copies of the ASEP system having configurations 770 and 771 which agree 
except at one site X, at which r]o(X) = and i]i(X) = 1, i.e., system has a hole and 
system 1 a particle. Allow this pair of systems to evolve under a coupled dynamics, so that 
attempts to jump from a given site to an adjacent one occur simultaneously. Then each 
time a jump occurs in either system the same jump occurs in the other, if possible; this 
synchronization can fail only when the extra particle in system 1 jumps, or when a particle 
in system jumps on the extra hole, and in these cases the mismatch position X will move. 
From a configuration (770,771) of this doubled system we may obtain a configuration r of 
the single ASEP with second class particle by taking r(X) = 2 and r(i) = rjo(i) = 771(1) 
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when % 7^ X; the dynamics for the doubled ASEP system corresponds to that described 
above of the system with a second class particle. Conversely, from a configuration r we 
may obtain two distinct ASEP configurations rjo and 771 by restricting attention to one or 
the other of the paired ASEP systems or equivalently by replacing the second class particle 
by respectively a hole or a first class particle; we will write 770 = ^o( T ) an d t]i = ^i(r). 
Similarly, from a measure A for the ASEP with second class particle we obtain ASEP 
measures \I/o(A) and ( A) giving the distribution of r) and rji under A. Clearly if X t is 
evolving under the second class particle dynamics then both \I/o(At) and \I/i(At) evolve 
under the simple ASEP dynamics. 

A measure A describing this system from the viewpoint of the second class particle, 
i.e., in a reference frame in which X = 0, was constructed explicitly in [17]; the construction 
is summarized in Section 4. This measure has spatial asymptotics corresponding to the 
shock, 

lim T~ k \ = v p± , (2.1) 

k— >±oo 

and is invariant under the natural dynamics for the system seen from the second class 
particle, under which the second class particle is always at the origin and a jump of this 
particle in the original dynamics becomes a jump of the rest of the system in the opposite 
direction. 

To obtain a measure on the original ASEP configurations from the measure A we may 
consider either \I/o(A) or \I/i(A). To be definite, let us focus for the moment on the former, 
which gives the distributions of the configuration rj , and denote it by jl; jl is obtained 
from A by replacing the second class particle at the origin with a hole. This is the measure 
referred to in the introduction. We may allow it to evolve to /it under the ASEP dynamics, 
or equivalently write fit — ^o(At); then viewed from the position X t it is time invariant: 

T~ x ^ t =% (2.2) 
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Thus ju furnishes an invariant description of the shock itself, ignoring its location. It follows 
from (2.1) that has the spatial asymptotics (1.3). 

Other points of view are possible. Thus one may obtain new descriptions of the shock 
by rather trivial shifts of viewpoint; for example, a constant one, to the second site to the 
right of the second class particle, a configuration-dependent shift, to the third empty site 
to its right, or a shift with additional randomness, such as a choice, with equal weights, 
between the two previous possibilities. 

An alternate measure for the shock is implicit in the description of the system with 
second class particle as a coupled pair of ASEP systems: the measure Jit = *i (At), obtained 
from At by replacing the second class particle by a first class particle, evolves with the ASEP 
dynamics, and is invariant in the sense of (2.2): T~ Xt Jlt = Jiq. We write ji = Jio = ^i(A); 
Ji is obtained from simply by replacing the empty site at the origin by a particle. Despite 
the fact that Jl and fl are equally valid as candidates for the description of the ASEP shock 
(in the terminology introduced in the next section, they are both invariant shock measures 
for the ASEP), it is not at all clear that they are shift equivalent, that is, differ by the sort 
of random change of viewpoint that we have been considering here. In Section 5 we will 
use the ideas of the next section to show that this is the case. 

Other choices may for some purposes be more tractable. The existence of a shock 
measure was first proved by Ferrari, Kipnis, and Saada [14] (FKS) using a random view- 
point Z t which Ferrari [15] later showed was related to X t by a random translation of finite 
mean. To construct Z t , consider again two copies of the ASEP system with configurations 
Co and Ci satisfying (o(k) < Ci(^) f° r au k, so that when there is a particle at site k in 
configuration (q there is also a particle at that site in configuration £i . Allow the system 
to evolve under the coupled dynamics described above, so that again sites k at which 
Co(k) — (i(k) = 1 and those at which (o(k) = and Ci(^) = 1 obey the dynamics of first 
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and second class particles, respectively. Let A* be a translation- and time-invariant mea- 
sure for this system in which the densities are given by (Co(0))a* = p~ and (Ci(0))a* = P+ 
(the existence of such a A* is established in [14]). Randomly select some second class par- 
ticle, i.e., some discrepancy between (o and (i (more precisely, condition on the presence 
of such a particle at the origin at time 0), and let Z t be its position at time t. At time 
t = 0, define an ASEP configuration r\ as follows: first, if Co(^) = Ci(^) = 1 then r](k) = 1, 
and if Co(^) = Ci(^) = then t](k) = 0; second, if (o(k) = and Ci(^) = 1) an d k is the 
j th site at which such a discrepancy occurs, counting from j = at Zq, then rj(k) = 1 with 
probability 1/(1 + (q/pY), and rj(k) = with the complementary probability, and all of 
these choices are independent. Allow the configuration rj to evolve with ASEP dynamics 
coupled to that of Co and Ci, so that if a jump occurs in any of the three systems it also 
occurs in any others in which it is possible. Then the distribution of rj at time t, viewed 
from the position Z t , is time independent and has the shock asymptotics (1.3). 

In all the examples considered so far the position of the viewpoint is not determined 
by knowledge of the ASEP configuration r)\ knowing r] gives only the distribution of the 
random viewpoint. This is an additional randomness beyond the configuration dependence 
discussed for the d = 2 Ising model in the introduction. Our last example is a construction 
of a viewpoint I on a non-TI measure p which, as in the d = 2 equilibrium example, is 
intrinsic. This means, first, that the viewpoint £(rj) depends only on the configuration r], 
with no additional randomness, and second, that the viewpoint behaves covariantly under 
translations, so that 

£(T V ) = £( V ) + 1. (2.3) 

The function £{rf) may be thought of as picking out a shock location in the configuration 
rj. Its intrinsic nature means that if v is any measure obtained from p by a shift of 
viewpoint then v and p look the same from the configuration dependent viewpoint £. In 

10 



the construction of £, \i need not be related to the ASEP dynamics; we require only that 
\x be a non-TI measure on ASEP configurations which has well-defined and configuration- 
independent asymptotic densities satisfying p + > p_: 



1 



N 



p± = lim 



TV^ioo \N\ + 1 



(2.4) 



k=0 



for ^-almost every r\. 

To define £{r\) we first define a function h v (j) to be the signed cumulative occupation 
from the origin to site j: 



The graph of h v has, for typical rj, slope p + (on a large scale) far to the right of the origin 
and p- far to the left. Now fix an irrational number satisfying p_ < p* < p + , and 
for any rj define £{rj) to be the integer j which minimizes h v (j) — p*j, whenever such a 
(necessarily unique) minimizing integer exists; for other rj, £(rj) is undefined. The special 
role played by the origin in the definition (2.5) of h v does not affect the value of £{rj). The 
construction is shown graphically in Figure 1. It is intuitively clear, and can be proved, 
that i is well defined and finite with probability one relative to p. The definition of I 
depends strongly on the choice of p*; slight changes in p* will cause large changes in the 
viewpoint £(rj) for some configurations rj. 

Note that the sets Sk = { V \ £{rf) = k } form a partition of the configuration space 
(up to a set of /(/-measure zero) which is nicely mapped by translations: T(Sk) = Sfc+i- 
Such partitions have been constructed for more general T in the context of discrete time 
dynamical systems by Gurevic and Oseledec [19]. 




if j>0, 
if j = 0, 
rj(i), if j < 0. 



(2.5) 
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The shock may look quite different from different viewpoints, in particular, from 
viewpoints I defined with different values of p*. In Figure 2 we show shock profiles (mean 
values (r)(k)) at site k relative to the viewpoint adopted) for the shock in the totally 
asymmetric (p = 1) model with densities p + = 0.7, p- = 0.2, seen from three different 
viewpoints: from the second class particle (that is, in the measure p) and from t\ and 
£2, the viewpoints constructed as above with p*i = %~ 1 p + + (1 — 7r _1 )p_ and p*2 = 
(l- 7 r- 1 )p + + 7r- 1 p- 

3. Equivalence of measures under random shifts 

In the previous section we have described implicitly an equivalence relation on prob- 
ability measures on the ASEP configuration space S = {0, 1} Z , under which two such 
measures p\ and P2 are equivalent if they differ by a configuration dependent random shift 
of viewpoint. Perhaps the simplest way to make this precise is in terms of a coupling for the 
two measures, that is, a measure on S x S with marginals p\ and p 2 on the first and second 
components. We say that p± and P2 are shift equivalent if there exists such a coupling p* 
and an integer- valued function Y such that for (771, 772) £ S x S, with 
/^-probability one; in this case we write 

p 2 = T~ Y p 1 . (3.1) 

Generalizations to measures on sets other than S on which a translation operator acts can 
easily be made. The coupling p* is sometimes referred to in the mathematical literature 
as a shift coupling [20,21] or an orbit coupling [22]. 

It is sometimes convenient to work with alternate formulations of this equivalence 
relation. One such is based on a variation of the well known Vasershtein distance [23] 
between two measures. Define ^(771 , 772) on S x S to be the minimum, over k, of the 
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number of sites at which 771 and T k rj 2 differ, and for measures \i\ and fi2 on S define 
D(ji U H2) by 

D(fjL 1 ,n 2 ) = w£(d} v ., (3.2) 

V* 

where the infimum is over couplings v* for fi\ and fi 2 - Then it can be shown [18] that fi\ 
and /j,2 are shift equivalent if and only if D(/j,i, fj, 2 ) = 0. The minimum in (3.2) is in fact 
achieved when v* = fi*, with fx* the coupling used to define (3.1). 

A second reformulation of the relation of shift equivalence is obtained by noting that 
if (3.1) holds then clearly 

fii(A) = fi2(A) for all A C S and T( A) = A, (3.3) 

that is, if A is translation invariant. Conversely, it can be shown [21,22] that if (3.3) holds 
then there exists a random position Y so that (3.1) is satisfied. For measures describing 
shocks, "interesting" TI sets A with nontrivial probability describe intrinsic properties 
of the shock; for example, we might take A to be the set of configurations 77 with a 
particle at the site three sites ahead of the position i(rj) defined in Section 2, so that 
fx{A) = ( V (e(r,) + 3) )„. 

Finally, as remarked earlier, when an intrinsic viewpoint £(rj) can be defined, that is, 
a function £(rj) defined almost everywhere with respect to both fi\ and 112 and satisfying 
(2.3), then ii\ and 112 are shift equivalent if and only if T~ i fii = T~ l fi 2 . 

The concept of shift equivalence allows us to give precise definitions of two natural 
concepts for shock measures in the ASEP (or similar systems). We say that a measure \i 
is an invariant shock measure if it has the spatial asymptotics (1.3) for some p± and if, 
when fi t is the measure evolving under the ASEP dynamics which satisfies fi$ = fx, fit is 
shift equivalent to fi for all t; for example, the measures ft and Ji of the previous section are 
invariant shock measures in this sense. We say that the ASEP has a unique shock measure 
(for given p±) if any two such invariant measures are shift equivalent. It seems natural to 
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conjecture that the ASEP has a unique shock measure in this sense for all p± satisfying 
p + > p_, but this has not been established. 

Let us now restrict attention to measures on S which, like the ASEP shock measures 
considered in Section 2, converge under spatial translation to distinct TI states. Suppose 
then that p± are translation invariant probability measures on S with p + ^ p_, and define 
a ramp measures to be a probability measure p on S which is asymptotic to p + to the 
right of the origin and to p- to the left: 

lim T~ k p = p±. (3.4) 

Associated to each ramp measures is a family of translation sums. These sums, under 
rather mild additional technical conditions on the measures involved [18], are invariant 
under a shift of viewpoint and furnish a complete characterization of shift equivalence. In 
a sense, the equivalence of (3.1) and (3.3) also provides a set of invariant quantities which 
determine the equivalence class of a measure p: the values p(A) for all TI sets A. The 
example of the translation invariant set given above, however, suggests correctly that these 
quantities are difficult to calculate. We will see in the next section that the translation 
sums are calculable for the ASEP shock measure. 

Suppose that p is a ramp measure and that / is a function on S which depends on only 
finitely many occupation numbers and which satisfies p+(f) = P-(f) = 0; for example, 
f(v) = ?7(1) — ^(0) is such a function and, if the asymptotic states are product measures, 
i.e., if p± = v p±1 then so is f{rj) = rj(k)(r](l) — rj(0)) whenever fc ^ 0,1. Then we may 
define the translation sum 

00 

M/) = E < t? V>m; ( 3 - 5 ) 

n= — oo 

A /x (/) will be finite for any ramp measure p for which the asymptotic behavior (3.4) is 
achieved with sufficient rapidity to guarantee that this sum converges. For example, if 
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f(rj) = rj(l) — i](0) then A M (/) = p + — p_ since (3.5) telescopes, and similarly A M (/) = 
(fl r )ju+ ~ ^ / = Tg — g for some #(77). However, we see no easy way to compute A M (/) 

for general // and /. 

The values of the translation sums characterize the (shift) equivalence classes of ramp 
measures in the following sense: under additional conditions describing the convergence at 
±00, two ramp measures p,\ and ^2 are shift equivalent if and only if 

A m (/) = A M2 (/) (3.6) 

for all / satisfying p+(f) = p_(f) = [18]. 

As a simple application of this result, consider a product measure v\ on S with density 
Pi(k) = (v(k)) Ul satisfying lim^-i-oo p\ (k) = p±, with p + 7^ p_; v\ is a ramp measure 
(asymptotic to u p+ and v p _) and the additional technical considerations needed for the 
above result are satisfied if the asymptotic limit is achieved sufficiently rapidly. Let z/ 2 be 
another such measure obtained by altering the density at the origin only: /O 2 (0) 7^ Pi(0), 
P2(k) = Pi(k) if k 7^ 0. Then v\ and v<i will typically have different translation sums and 
hence not be shift equivalent; for example, consideration of the translations sums for the 
functions fk(v) = v(k)(v(l) ~ ^(0)), k > 1, shows that v\ and z/ 2 will not be shift equivalent 
unless pi(k) + p\{—k) is independent of k. 
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4. Translation sums for the ASEP shock measure 

In this section we show how to compute the translation sums A~(/) for the ASEP 
shock measure ju which, as discussed in Section 2, is obtained from the invariant shock 
measure A for the system with second class particle at the origin by replacing that particle 
by a hole. In this case the asymptotic measures fx± are the product measures v p± . We will 
continue to denote a typical ASEP configuration by 77 (rj(i) = or 1) and a configuration 
in which second class particles may occur by r (r(i) = 0, 1, or 2); in this section such a r 
will always contain a single second class particle located at the origin. 

In [17] (which was an extension to the general ASEP of the results of [16] for the 
totally asymmetric model, in which p = 1) it was shown that the measure A can be written 
in terms of two vectors \v) and (w\ and three operators A, D and E satisfying the following 
algebraic rules: 



pDE - qED = 


(P- 


«)[(!- P-)(l 


- p+)D + p-p+E], 


(4.1) 


pAE - qEA = 


(p- 


«)(1-P-)(1- 


-p+)A 


(4.2) 


pDA - qAD = 


(P- 


q)p+P-A, 




(4.3) 


(D + E)\v} = 


\v), 






(4.4) 


(w\(D + E) = 


H; 






(4.5) 


(w\A\v) = 


1. 






(4.6) 



Specifically, the probability of the set of configurations specified by the occupation numbers 
C(z) (C(z) = 0, 1) of m consecutive sites to the left of the second class particle (which is 
located at the origin) and n consecutive sites to its right, 

A({ r I r(i) = C(i). i = -m, . . . , -1 and i = 1, . . . , w }) = P m ,„(C), (4.7) 

can be written as the matrix element in which a first class particle is represented by a 
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matrix D, an empty site by a matrix E, and the second class particle by a matrix A: 

Pm,n(0 = w\ J} [C(^+(l-CW)^U|n[C(j^+(l-CW)^]|k)- (4-8) 
Ki=— m ) {j^ 1 J 

For example, the probability of finding occupation numbers 10 1 immediately to the left 
and 110 immediately to the right of the second class particle, that is, of the local 
configuration 10120110 0, is given by (w\DEDAED 2 E 2 \v). 

It follows from the above that the microscopic shock profile, defined as the average 
occupation (t(to))~ = {r]{n))^ at site n ^ 0, is given for n > by 

(r(n))~ = (w\A(D + E^D^). (4.9) 

The exact expression of this profile was given in equations (4.1-4.5) of [17], where it was 
also shown that the profile has the symmetry property 

(r(n)) 3: +(T(-n))~ = p + +p-. (4.10) 

Let us call a finite product of D's and E : s a word. To every function of k consecutive 
site occupation numbers, say f(rj(l),rj(2),...,rj(k)), there is naturally associated a linear 
combination of words of length k; for example, if / = ??(1)(1 — ??(3)) + 377(5)77(6) then 

W = D(D + E)E(D + Ef + 3(D + E) 4 D 2 . (4.11) 

We denote the word of length zero by 1 and associate it with the function / = 1. For 
every linear combination W of words we will define below a number r(VF) and will then 
develop, from the algebraic rules (4.1)-(4.6), new rules which enable us to calculate r(VF). 
This will determine the translation sums A^(/) defined in (3.5), since we will show below 
that when / and W are related as above and in addition (/) M± = 0, A^(/) = r(VF). 
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Suppose then that W is a linear combination of words and define 



F(W)= lim {T L M (W) - Lr+(W) - Mr_(W)}. (4.12) 



Here 



T l ,m(W)= ^(w\(D + E) l W(D + E) m \v) 



(4.13) 



9=0 



with 



D = D, E = E + OA, (4.14) 



and W the operator obtained from W by replacing each D and E in W by D and E; 
r±(W) denotes the number obtained by replacing each D and E in W by p± and 1 — p± 
respectively. It is clear from (4.13) and the definition (4.14) of D and E that the right 
hand side of (4.13) is a sum of matrix elements of products of the operators D, E, and A, 
with each product containing a single operator A; each such matrix element is computable 
from (4.1)-(4.6) and represents a probability calculated in the measure A. Moreover, the 
limit in (4.12) exists because A converges exponentially fast to v p± at ±oo [17]. Thus T(W) 
is well defined. 

It follows directly from the definition of T(W) that 

T{a 1 W l + a 2 W 2 ) = a^Wt) + a 2 r(W 2 ), for any a x and a 2 ; (4.15) 

T(W[D + E]) = T(W)+r_(W); (4.16) 

T([D + E]W) = T(W) +r + (W); (4.17) 

r(l) = 0. (4.18) 

Also, one can easily show that D and E satisfy (4.1), and hence for any W\ and W 2l 



T(W 1 (pDE - qED)W 2 ) = (p - g)r(Wi[(l - p_)(l - p+)D + p_p+E]W 2 ). (4.19) 
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The five rules (4.15)-(4.19), together with the value of T(D), which we will derive 
below, allow one to calculate r(W) for any linear combination of words W. The argument, 
based on a recursion on the length of a word W, is almost identical to that given in [24,17] 
to show that (4.1)-(4.6) are sufficient to calculate any matrix element of the form (4.8). 
Consider, for example, T(DW n ), where W n is a word of length n > 1 with k factors E: 

p n T(DW n ) = q k p n - k T(W n D) + l.o.t. 

= —q k p n ~ k T(W n E) + q k p n - k r-(W n ) + l.o.t. 
= -q n T(EW n ) + q k p n - k r-(W n ) + l.o.t. 

= q n T(DW n ) - q n r + (W n ) + q k p n - k r_(W n ) + l.o.t. (4.20) 

Here l.o.t. (lower order terms) denotes linear combinations of r( W) for words W of length 
n or less. Since q < p, (4.20) can be solved for T(DW n ). For example, since (4.16) and 
(4.18) imply that F(D + E) = 1, 

pT{D 2 ) = - P T(DE) +pT(D) +pp_ 

= -qT(ED) - (p - g)r((l - p+)(l - p-)D + P-P+E) +pT(D) + pp_ 

= qT(D 2 ) + pp. - qp+ - (p - q)p_p + + (p - q)(p+ + P -)T(D), (4.21) 

and hence 

r(D 2 ) = PP ~ - QP+ - P+P- + ( P+ + P-)T(D). (4.22) 
p-q 

If W is a word of length k which contains j factors E and if for i = 1, . . . ,j, 
is the operator product obtained from W by replacing the i th factor E by A, then from 
(4.4), (4.5), (4.13), and (4.14), 

L-l j M-1 

F LjM (W) = ^2(w\A(D + EyW\v) + J2(w\W {i) \v) + {w\W(D + E) l A\v). (4.23) 

i=0 i=l i=0 
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More generally, if f(rj(l), 77(2) , rj(k)) is a function of the occupation numbers of sites 
l,...,k and W the corresponding linear combination of words of length k (see (4.11)), 
then (4.23) implies that 



L-l 



r L , M (w)= £ (T7> A ; (4.24) 



i=-M-k 



Note that, as is clear in (4.23), the averages are taken in the measure /i obtained from A 
by replacing the second class particle at the origin by a hole, that is, taking r](0) = 0. For 
example, if / = 7y(l)(l - 77(3)) + 377(5)77(6) then from (4.11) and (4.24), 

L 

r L , M (D(D+E)E(D+E) 3 + 3(D+E) 4 D 2 ) = Mi)(l-»7(i+2))+377(i+4)77(i+5))£. 

(4.25) 

Equations (4.24) and (4.10) imply that 

F(D) = ^lim I M% - L (P+ + P-) J = (^(0)>- = 0. (4.26) 

Moreover, (4.24) leads immediately to the calculation of translation sums A^-(/), for if / 
is a function satisfying (/) M± = and is the corresponding linear combination of words, 
then (4.12), (4.24), and the relations (f)„ = r±(W) imply that 

A~(f) = T(W). (4.27) 

The translation sums are in some cases related to expectation values in the measure 
/i in a surprisingly simple way. For example, let fi(rj) be a function of rj(—m), . . . , rj(—l) 
and f2(v) a function of 77(2), . . . ,?7(n), and let h(rf) = 77(1) — 77(0). Then n±(fihf 2 ) = 0, 
so that A^-(/i/t,/ 2 ) is defined. We will show that 



A~(M/ 2 ) = (P+ - P-){fi(T- 1 f2))~. (4.28) 
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The operator corresponding to h is A' = (D + E)D - D(D + E) = ED - DE; in view 
of (4.8) and (4.27), (4.28) will follow if we show that for any W u W 2 , 

T(WiA'W 2 ) = (p+ - p-)(w\W 1 AW 2 \v). (4.29) 

To verify (4.29) we show that, up to a normalization, I^WiA'T^) is determined by the 
same algebraic rules (4.1)-(4.6) which determined (u^Wi-AW^u)- First, D and E of course 
satisfy (4.1), and this in turn implies that A' satisfies the analogue of (4.2) and (4.3): 

pA'E - qEA' = (p - q)(l - p_)(l - p+)A', (4.30) 
pDA' - qA'D = (p - q)p+p-A'. (4.31) 

Moreover, since r + (W 1 A'W 2 ) = r_(WiA'W 2 ) = for any W 1 ,W 2 , (4.16) and (4.17) give 

T([D + E]W 1 A'W 2 ) = YiWxA'W^D + E]) = T(WiA / W 2 ), (4.32) 

which is the analogue of (4.4)-(4.5). Finally, the analogue of (4.6) is 

F(A') = F([D + E]D) - F(D[D + E}) = p+ - p- (4.33) 

Since the rules (4.1)-(4.6) determine all matrix elements (tu|W r i^4W"2|'i;), (4.1) and (4.30)- 
(4.33) imply that these will agree with the r^iA'Vl^) up to a normalization determined 
by comparing (4.6) and (4.33), that is, (4.29) holds. 
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5. Consequences of the calculation of the translation sums 



In this section we apply the calculation of the translation sums A~(/) outlined in 
Section 4 to discuss the shift equivalence of various measures describing the ASEP shock 
and to identify certain intrinsic features of the ASEP shock. 

In Section 2 we observed that the measure Jl, obtained from A by replacing the second 
class particle at the origin by a first class particle (or equivalently from fi, by replacing the 
hole at the origin by a particle) , stood on an equal footing with in providing a description 
of the shock. We now show that 

A~(/) = A~(J) (5.1) 

for any / satisfying (/) /i± = 0, so that by the general results referred to in Section 3, ju 
and fx differ only by a shift of viewpoint (which in general will be a random shift with 
distribution depending on the configuration). 

For suppose that, instead of using (4.14), we had defined D and E by D = D + 6 A 
and E = E, arriving at new quantities T(W). Then all the considerations of Section 4 
would have been essentially unchanged, except that averages in (4.24) would be computed 
in the measure Ji; as a result, (4.26) would become f(D) = 1. This would lead to 

r+(W)-r-(W) 



T(W) = T(W) + [T(W) - r(w)]- 



p+ - p- 



+ (82) 

P+ - P- 



this equation may be verified by noting that ri(W) = T(W) - T(W) and T 2 {W) = 
r + (W) — r-(W) both satisfy (4.15), (4.18), (4.19), and the homogeneous versions of (4.16) 
and (4.17), and that (5.2) holds for W = D, so that the reduction procedure (4.20) 
yields (5.2) for all W. As a consequence, from (4.27) and the corresponding equation 
A~(/) = r(W), the translation sums satisfy (5.1). 
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We next turn to the calculation of the translation sums associated with the particular 
family of functions f n {rj) = (77(0) — p_) (rj(n) — p + ), n > 1. Let us write 

00 

$ n = A~(f n ) = <fa(0 ~ P-)(v(i + n) - (5-3) 

i= — 00 

Since / n (7/) - / n+1 (7/) = (t/(ti + 1) - 7/(n))(p_ - 77(0)) = (T n F)(rj), where Ffa) = (77(1) - 
//(0))(p_ - rj(-n)), (4.28) implies that 



$ n - $ n+1 = ( P+ - p_)[p_ - (v(-n))~] = (p+ - p-)[(ri(n))~-p + }, (5.4) 



where we have used the symmetry (4.10) of the profile. Thus 



$ n - lim = (p+ - p_) Y^[(v(i)h- P+l- (5-5) 



But $1 can be evaluated from (4.15)-(4.18), (4.22), and (4.26): 



$1 = (p+ - p-) x>(% -p+}- = pp ~ qp+ - p+p- (5.6) 



The right hand side of (5.5) in the case n = 1 can be evaluated from formulas (4.1)-(4.5) 
of [17] and shown also to be given by (5.6); thus limfc^oo $fc = (as one would expect) 
and hence 



n = (p+-P-)£[(r/(%-p + ]. (5.7) 



To understand the consequences of equation (5.7) we recall [17] that in general the 
profile (r}(n))^ decays exponentially to its asymptotic value: for n ^> 1, 



fa(n)) ; -p + ~Cn^- om , (5.8) 
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where C > if q/p < x*, C = (and in fact (rj(n))~ = p + for all n > 0) if g/p = x*, and 
C < if q/p > x*, with x* = — p+)/p+(l — p_) (explicit values of a and 7 are given 
in [17]). Thus (5.7) implies that, except on the line q/p = x*, $ n decreases exponentially 
to as n — > 00, with characteristic length a -1 . Since the $ n are intrinsic properties, this 
gives an intrinsic characteristic size for the shock. 

A somewhat surprising aspect of the shock profile, found in [16] for p = 1 and in [17] 
for all p satisfying q/p < x*, is the "overshoot": (r)(n))^ > p + for n > 1 (corresponding 
to C > in (5.8)). A natural question is whether this overshoot is a by-product of the 
choice to view the shock from the second class particle, and might be eliminated by the 
adoption of another viewpoint. While we cannot answer this completely, (5.7) implies that 
for q/p < x* there is no viewpoint from which the shock measure is described by a product 
measure with density increasing monotonically from p_ to p+, since such a measure would 
lead to a negative <£> n . Similarly, for q/p 7^ x* there is no point of view from which the 
shock measure would be a product measure with density p+ to the right of the origin and 
density p_ to the left of the origin, since for such a measure, all the $ n would vanish. 

6. Shocks in other ASEP models 

Several different models, with ASEP dynamics but with specific initial conditions, 
boundary conditions or minor modifications, have been shown to give rise to shocks: 

1. An infinite system with, as initial condition, a product measure with density p+ 
to the right and p- to the left of the origin. This is the case discussed in the 
introduction and in [14,15,17]. 

2. A system with ring geometry and a blockage bond, at which the jump rate is less 
than that at other bonds in the system [25]. 

3. A system with ring geometry and a slow second class particle [26,27]. 
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4. A system with open boundary conditions on its first order transition line, discussed 
below and in [28]. 

Note that in examples 2 and 3 the shock in the ring geometry is caused by the blocking 
bond or particle but is located far from it, in a region where the usual ASEP dynamics 
holds. We believe that, when the size of the system goes to infinity, all of these different 
models lead in fact to shocks which are shift equivalent, i.e., that one may choose for each 
shock a viewpoint such that, seen from this viewpoint, all the shocks are described by the 
same ramp measure. (The relation of shocks observed in similar but more complicated 
systems [29] to those of the ASEP remains to be investigated.) This equivalence is not 
obvious a priori but can be verified, at least in cases where exact expressions permit exact 
computations, by showing that all the translation sums A(/) are identical to those found 
in Section 4. 

We illustrate this by computing the translation sums for the ASEP with open bound- 
ary conditions, in the totally asymmetric (p = 1) case. This is a system with N sites, in 
which particles enter the system at site 1 with rate a and escape at site N with rate f3. 
In [28], it was shown that the steady state of this system can be fully described by an 
algebra rather similar to the one discussed in Section 4 (the weight of a configuration can 
be written as the matrix element of a matrix product where a matrix D is used when a site 
is occupied and a matrix E when the site is empty). Using this algebra, in principle, any 
expectation can be calculated for systems of any size and for any values of the parameters 
a and (3. In particular, it was shown that if a < min{/3, 1/2} the system is in a low density 
phase, described in the bulk by a product measure with density a, while if (3 < min{a, 1/2} 
the system is in a similar high density phase, with bulk density 1 — (3. 

On the line a = j3 < 1/2 there is a first order phase transition between these two 
phases. When the system is on this phase coexistence line one can show [28] that for large 
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N and for % far from the boundaries (j>l and A^-i > 1) the profile is linear in i/N, 



and by a similar computation, not given in [28] , that nearest neighbor two point correlations 
behave in the same way: 



The interpretation of this behavior is that there is a shock in the system, separating a 
region of low density a to its left and of high density 1 — a to its right, and that the 
linear dependence on i/N arises because the shock is equally likely to be anywhere in the 
system (apart from some corrections of order 1/N due to boundaries). It seems natural 
to suppose as well that the structure of the shock is independent of its location, so that if 
we determine the shock position in the configuration rj by some function £(rj) as described 
in Section 2, and write \i m for the system measure conditioned on l{yj) = m, then the 
measures fx m are just translates of one another, \i m = T m ~ m ' fi m ', for N large and m far 
from the boundaries. 

If such an interpretation is correct, then the expectation of any function of the occu- 
pation numbers near site i should for large N show linear behavior in i/N similar to that 
of (6.1) and (6.2). The leading terms give no information about the shock structure, but 
it follows also from the picture described above that there is a probability of order 1/N 
that the shock is located near site i, so that to describe the shock properties one must 
compute the terms of order 1/N in expressions like (6.1) and (6.2). However, these order 
1/N terms may contain not only contributions coming from the events where the shock is 
in the neighborhood of i, but also contributions from boundary effects, in particular from 
the precise definition of the coordinate i. 




(6.1) 




(6.2) 
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We now show that we may calculate these order 1/N terms and thus obtain additional 
support for this picture. If / = f(r](i + 1), r](i + 2), ...rj(i + k)) is a function depending on 
the occupation numbers of sites between % + 1 and % + k, with k <C N, then it can be shown 
that for i far enough from the boundaries the expectation of / takes the form 

(/) = WTb + NTb + {N) ■ (6 - 3) 

Here b is a constant of order 1, independent of / and i, and W is the linear combination 
of words of length k associated to /, as in Section 4 (see (4.11)). Equation (6.3) defines 
7(W); the correct choice of 6, which may be regarded as representing the boundary effects, 
guarantees that 7(W) is independent of i. Equation (6.3) is trivial for / = 1 (corresponding 
to k = 0) and leads to 

7 (1) = 0. (6.4) 

For / = rj(i), (6.3) may be verified from the asymptotics derived in [28]; this calculation 

leads to explicit values for b and of 7(-D), but the values of these parameters are not needed 

for the calculation of the translation sums. For other words W one can use the algebra of 

[28] to verify (6.3) by induction on the length of W, and also show that the 7(W) have 

the following properties: 

l(W(D + E)) = 1 (W) + r-(W) (6.5) 

<y((D + E)W) =<y(W) +r + (W) (6.6) 

7(aiWi + a 2 W 2 ) = ai7(Wi) + a 2 -f(W 2 ) (6.7) 

^(WrDEW^ = a(l - a)[ 1 {W l DW 2 ) + ^{WxEW^} (6.8) 

Here r + and r_ are defined with the densities p + = 1 — a and p_ = a. Thus the 7(W) 
satisfy the same rules as r(W) of Section 4, and one can conclude (see the argument 
following (5.2)) that the general expression of these 7(W) is given by 

7 W = r W + 7 (D) ^ +W ~ /A - W (6.9) 

p+ - p- 
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Now suppose that the function / satisfies (f) Va = (f) Ul _ a = 0, so that (6.3) becomes 

(f) = J ^p- + o(l/N). (6.10) 

This result has a simple interpretation in terms of the heuristic picture, described above, 
that the steady state of the open system contains a shock equally likely to be located at any 
site and with structure independent of its location. Thus in computing (/) we are averaging 
the expected value in the measure conditioned on the shock position being m, (/) Mm , over 
all N possible shock positions. Since the measures \i m differ only by translation this is, up 
to a factor 1/iV, just the calculation of the translation sum A Mm (/) for any fixed m (N 
must be large enough so that \i m is effectively an infinite volume measure). Thus (6.10) 
implies that A Mm (/) = 7(W). But (6.9) implies that for such a function /, 7(W) = r(W), 
so that from (4.27), A Mm (/) = A~(/): the translation sums for the shock in the open 
system are the same as those calculated in Section 4. Thus these two shocks are shift 
equivalent. 
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Figure captions 



Figure 1. Construction of the intrinsic shock location function £(rj). 

Figure 2. Shock profiles for p = 1, p + = 0.7, p- = 0.2, from three different viewpoints: 
(i) the second class particle (diamonds, solid line), (ii) t defined with p* = ir~ 1 p + + (1 — 
n~ 1 )p- (plusses, dashed line), (iii) £ defined with p* = (1 — ir~ 1 )p + + n~ 1 p- (squares, 
dotted line). 
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